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At LHC energies, perturbative QCD (pQCD) provides a powerful calculational tool. For Z and 
W transverse momentum spectra, pQCD theory agrees with the CDF[[j]] and DO |Q] data very well at 
Tevatron energies^]. The LHC pp program will test pQCD at an unprecedented energy. The heavy-ion 
pq" ■ program at the LHC will make it possible for the first time to observe the full p t spectra of heavy vector 

bosons in nuclear collisions and will provide a testing ground for pQCD resummation theory [ffl]. 

O ■ 

In nuclear collisions, the power corrections will be enhanced by initial and final state multiple 
q | scattering. As we will show, the high-twist effects are small at LHC for heavy boson production. The 

only important nuclear effect is the nuclear modification of the parton distribution function (shadowing). 
Because of their large masses, W and Z° will tell us about the nuclear Parton Distribution Function 
(nPDF) at large scales. Since for the p t spectra of heavy bosons contributions from different scales need 
to be resummed, the p t spectra for heavy vector bosons will provide information about the evolution 
of nPDFs from small scales to large scales. Since it is more difficult to detect W than Z°, we will 
concentrate on discussing Z° here (the results for production are very similar[B|]). 

m ■ 

Resummation of the large logarithms in QCD can be carried out either in p y-space directly, or 
CN \ in the so-called "impact parameter", 6-space, which is a Fourier conjugate of the p y-space. Using the 

renormalization group equation technique, Collins, Soper, and Sterman (CSS) derived a formalism for 
the transverse momentum distribution of vector boson production in hadronic collisions^]. In the CSS 
formalism, non-perturbative input is needed for the large b region. The dependence of the pQCD results 
Q_i' on the non-perturbative input is not weak if the original extrapolation proposed by CSS is used. Recently, 

a new extrapolation scheme was introduced, based on solving the renormalization group equation includ- 
ing power corrections^]. Using the new extrapolation formula, the dependence of the pQCD results on 
the non-perturbative input was significantly reduced. 



For vector boson (V) production in a hadron collision, the CSS resummation formalism yieldsO 
"*-!■ da(h A + h B ^V + X) _ 1 [, 27 ci p T .t. 



dM 2 dydp 2 T (2vr)' 



J d 2 be l P Trb W(b,M,x A ,x B ) + Y(p T ,M,x A ,x B ) , (1) 



where x A = e y Mj yfs and x B = e~ y M/y/s, with rapidity y and collision energy y/s. In Eq. (|J), the W 
term dominates the p x distributions when pj- < M, and the Y term gives corrections that are negligible 
for small p t, but become important when p t ~ M. 

The function W(b, M, x A ,x B ) can be calculated perturbatively for small b, but an extrapolation to 
the large b region requiring nonperturbative input is necessary in order to complete the Fourier transform 
in Eq. ([[]). In oder to improve the situation, a new form was proposed[|3|] by solving the renormalization 
equation including power corrections. In the new formalism, W(b, M, x A ,x B ) = W pert (b, M, x A ,x B ), 
when b < b max , with 

W per \b, M,x A ,x B )= e s $M w(b, c/b, x A ,x B ) , (2) 

where all large logarithms from ln(l/6 2 ) to ln(M 2 ) have been completely resummed into the exponential 
factor S(b, M), and c is a constant of order unity [Q]. For b > b max , 

W(b,M,x A ,x B ) = W pert (b 

max) 1 \ u > u max ) , 0) 

where the nonperturbative function F NP is given by 

F np = exp{ _ l n{ M 2 b 2 max /c 2 ) [ 91 ((b 2 r - (blaxT) +92 (b 2 - b 2 max )] - g 2 (b 2 - b 2 max )}. (4) 



Here, b max is a parameter to separate the perturbatively calculated part from the non-perturbative input. 
Unlike in the original CSS formalism, W(b, M, xa,xb) is not altered, and is independent of the non- 
perturbative parameters when b < b max . In addition, the 6-dependence in Eq. (|j) is separated according 
to different physics origins. The (b) a -dependence mimics the summation of the perturbatively calcu- 
lable leading power contributions to the renormalization group equations to all orders in the running 
coupling constant a s ([i). The 6 2 -dependence of the g2 term is a direct consequence of dynamical power 
corrections to the renormalization group equations and has an explicit dependence on M. The g 2 term 
represents the effect of the non-vanishing intrinsic parton transverse momentum. 

A remarkable feature of the 6-space resummation formalism is that the resummed exponential 
factor exp[S(6, M)) suppresses the 6-integral when b is larger than 1/M. It can be shown using the 
saddle point method that, for a large enough M, QCD perturbation theory is valid even at p t = 0[Q]. As 
discussed in Ref.s [^, |5|, the value of the saddle point strongly depends on the collision energy y/s, in 
addition to its well-known M 2 dependence. The predictive power of the 6-space resummation formalism 
should be even better at the LHC than at the Tevatron. 

In Z° production, since final state interactions are negligible, power correction can arise only 
from initial state multiple scattering. Equations (|3|) and (Q) represent the most general form of W, and 
thus (apart from isospin and shadowing effects, which will be discussed later), the only way nuclear 
modifications associated with scale evolution enter the W term is through the coefficient g 2- 

The parameters g\ and a of Eq. (|j) are fixed by the requirement of continuity of the function W(b) 
and its derivative at b = b max . (The results are insensitive to changes of b max € [0.3 GeV -1 ,0.7 GeV -1 ]. 
We use bmax — 0-5 GeV ^.) The value of gi and c/2 can be obtained by fitting the low-energy Drell-Yan 
data. These data can be fitted with about equal precision if the values g 2 = 0.25 ± 0.05 GeV 2 and 
52 = 0.01 ± 0.005 GeV 2 are taken. As the b dependence of the 52 and <?2 terms in Eq. (^) is identical, 
it is convenient to combine these terms and define G2 = ln(M 2 b 2 nax /c 2 )g2 + 52 • Using the values of 
the parameters listed above, we get G2 = 0.33 ± 0.07 GeV 2 for Z° production in pp collisions. The 
parameter G2 can be considered the only free parameter in the non-perturbative input in Eq. @, arising 
from the power corrections in the renormalization group equations. An impression about the importance 
of power corrections can be obtained by comparing results with the above value of G 2 to those with 
power corrections turned off (G2 = 0). We therefore define the ratio 



The cross sections in the above equation and in the results presented here have been integrated over 
rapidity (—2.4 < y < 2.4) and invariant mass squared. For the PDFs, we use the CTEQ5M set[^]. 

Figure 1 displays the differential cross sections and the corresponding R g 2 rat i° (with the limiting 
values of G2 = 0.26 GeV 2 (dashed) and G2 = 0.40 GeV 2 (solid)) for Z° production as functions of p t 
at y/s = 14 TeV The deviation of Rg 2 from unity decreases rapidly as p t increases, and it is smaller 
than one percent for both y/s = 5.5 TeV (not shown) and y/s = 14 TeV in pp collisions, even when 
p t = 0. In other words, the effect of power corrections is very small at the LHC for the whole p t 
region. 

Without nuclear effects on the hard collision, the production of heavy vector bosons in nucleus- 
nucleus (AB) collisions should scale as the number of hard collisions, AB. However, there are several 
additional nuclear effects on the hard collision in a heavy-ion reaction. First of all, the isospin effect, 
which come from the difference between the neutron PDFs and the proton PDFs, is about 2% at LHC. 
This is expected, since at the LHC x ~ 0.02, where the u — d asymmetry is very small. 

The dynamical power corrections entering the parameter g 2 should be enhanced by the nuclear 
size, i.e. proportional to A 1 / 3 . Taking into account the A-dependence, we obtain G2 = 1.15 =t 0.35 
GeV 2 for Pb+Pb reactions. We find that with this larger value of G2, the effects of power corrections 
appear to be enhanced by a factor of about three from pp to Pb+Pb collisions at the LHC. Thus, even the 
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Fig. 1: (a) Cross section da/dpr for Z° production in pp collisions at the LHC with y/s = 14 TeV; (b) Rg 2 defined in Eq. (g) 
with G 2 = 0.26 GeV 2 (dashed) and 0.40 GeV 2 (solid). 

enhanced power corrections remain under 1 % when 3 GeV < px < 80 GeV. This small effect is taken 
into account in the following nuclear calculations. 

Next we turn to the phenomenon of shadowing, expected to be a function of x, the scale //, and of 
the position in the nucleus. The latter dependence means that in heavy-ion collisions, shadowing should 
be impact parameter (6) dependent. Here we concentrate on impact-parameter integrated results, where 
the effect of the 6-dependence of shadowing is relatively unimportant [Q], and we focus more attention on 
scale dependence. We therefore use EKS98 shadowing^] in this work. We define 



where Za and Zb are the atomic numbers and A and B are the mass numbers of the colliding nuclei, 
and the cross section da(px, Za/A, Zb / B)/dpx has been averaged over AB, while da(px)/dpT is the 
pp cross section. We have seen above that shadowing remains the only significant effect responsible for 
nuclear modifications. 

In Fig. 2(a), R s h (solid line) is surprising, because even at p j> = 90 GeV, x ~ 0.05, and we are 
still in the "strict shadowing" region. Therefore, the fact that R s h > 1 for 20 GeV < px < 70 GeV is 
not "anti-shadowing". To better understand the shape of the ratio as a function of p x, we also show R s h 
with the scale fixed at the values 5 GeV (dashed line) and 90 GeV (dotted), respectively, in Fig. 2(a). 
The nuclear modification to the PDFs is only a function of x and flavor in the calculations represented 
by the dashed and dotted lines. These two curves are similar in shape, but rather different from the solid 
line. In b space, W(b, M, xa,xb) is almost equally suppressed in the whole b region if the fixed scale 
shadowing is used. However, with scale-dependent shadowing, the suppression increases as b increases, 
as a result of the scale \i ~ 1/6 in the nPDF We can say that the scale dependence re-distributes the 
shadowing effect. In the present case, the re-distribution brings R s h above unity for 20 GeV < px < 70 
GeV. When px increases further, the contribution from the Y term starts to be important, and R s h dips 
back below one to match the fixed order pQCD result. 

We see from Fig. 2 that the shadowing effects in the px distribution of Z° bosons at the LHC are 
intimately related to the scale dependence of the nPDFs, on which we have only very limited data [|8[] . 
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Fig. 2: Cross section ratios for Z° production in Pb+Pb at y/s = 5.5 TeV: (a) R a h of Eq. (g) (solid line), and R s h with the 
scale fixed at 5 GeV (dashed) and 90 GeV (dotted); (b) Rg 2 of in Eq. (|) with G 2 = 0.8 GeV 2 (dashed) and 1.5 GeV 2 (solid). 

Theoretical studies (such as EKS98) are based on the assumption that the nPDFs differ from the parton 
distributions in the free proton, but obey the same DGLAP evolution^. Therefore, the transverse mo- 
mentum distribution of heavy bosons at the LHC in Pb+Pb collisions can provide a further test of our 
understanding of the nPDFs. 

In summary, higher-twist nuclear effects appear to be negligible in Z production at LHC en- 
ergies. We have demonstrated that the scale dependence of shadowing effects may lead to unexpected 
phenomenology of shadowing at these energies. Overall, the Z transverse momentum distributions can 
be used as a precision test for leading-twist pQCD in the TeV energy region for both, proton-proton and 
nuclear collisions. We propose that measurements of Z spectra be of very high priority at the LHC. 

Acknowledgments 

This work was supported in part by the U.S. Department of Energy under DE-FG02-86ER-40251. 
References 

[1] T. Affolder et al, CDF Collaboration, Phys. Rev. Lett. 84, 845 (2000). 



[2] B. Abbott et al, DO Collaboration, |hep-ex/0010026| (2000) 



[3] J.W. Qiu and X.F Zhang, Phys. Rev. Lett. 86, 2724 (2001); Phys. Rev. D 63, 114011 (2001). 

[4] J.C. Collins, D.E. Soper and G. Sterman, Nucl. Phys. B 250, 199 (1985). 

[5] X.F. Zhang and G. Fai, Phys. Rev. C 65, 064901 (2002). 

[6] H. L. Lai et al., Eur. Phys. J. C 12, 375 (2000). 

[7] Y. Zhang, G. Fai, G. Papp, G.G. Barnafoldi, and P. Levai, Phys. Rev. C 65, 034903 (2002). 

[8] KJ. Eskola, V.J. Kolhinen, and C.A. Salgado, Eur. Phys. J. C 9, 61 (1999). 



